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Small-angle, two-particle correlation functions have been 
measured for 36 Ar+ 112 ' 124 Sn collisions at E/A = 61 MeV. To- 
tal momentum gated neutron-proton (np) and proton-proton 
(pp) correlations are stronger for the 124 Sn-target. Some of 
the correlation functions for particle pairs involving deuterons 
or tritons (nd, pt, and nt) also show a dependence on the 
isospin of the emitting source. 

PACS number(s): 25.70.Pq 

1 Corresponding Author: Roberta Ghetti, Department 
of Physics, Lund University, Sweden E-mail address: 
roberta.ghetti@nuclear.lu.se 

The isospin dependence of the nuclear equation of 
state (EOS) is probably the most uncertain property of 
neutron-rich matter. This property is essential for the 
understanding of extremely asymmetric nuclei and nu- 
clear matter as it may occur in the r-process of nucle- 
osynthesis or in neutron stars [1]. In order to study the 
isospin-dependent EOS, heavy ion collisions with isotope 
separated beam and/or target nuclei can be utilized [2]. 
In these collisions, excited systems are created with vary- 
ing degree of proton-neutron asymmetry. A noticeable 
isospin dependence of the decay mechanism has been 
predicted [3-7]. Sensitive observables should be: pre- 
cquilibrium neutron/proton emission ratio [8], isospin 
"fractionation" [9-12], isoscaling in multifragmentation 
[13], neutron and proton flows [14]. 

Recently, the two-nucleon correlation function has 
been considered as a probe for the density dependence 
of the nuclear symmetry energy [15,16]. In this theo- 
retical study with an isospin-dependent transport model 
(IBUU), it was shown that a stiff EOS causes high mo- 
mentum neutrons and protons to be emitted almost si- 
multaneously, thereby leading to strong correlations. A 
soft EOS delays proton emission, which weakens the np 
correlation. In this paper we study experimental two- 
particle correlation functions of systems similar in size, 
but with different isospin. This work shows that, indeed, 
an isospin signal can be derived. 



Two-particle correlation functions were measured in 
E/A = 61 MeV 36 Ar-induced collisions on isotope- 
separated targets of 112 Sn and 124 Sn. The experiment 
was performed at the AGOR Superconducting Cyclotron 
of KVI (Groningen). The interferometer consisted of 16 
CsI(Tl) detectors for light charged particles, mounted at 
a distance 56-66 cm from the target in the angular range 
30°<#<114°, and 32 liquid scintillator neutron detectors, 
mounted 2.7 m from the target behind the "holes" of the 
Csl array, in matching positions to provide the np in- 
terferometer [17]. In this analysis, only data from the 
angular range 60 <6'<120 were used for neutrons. Fi- 
nally, 32 phoswich modules from the KVI Forward- Wall 
were mounted in the angular range 6°<#<18°, to collect 
information on the centrality of the collision. At least one 
fragment in the Forward- Wall was always required in our 
selected events, which biases our data towards midpe- 
riphcral collisions [18,19]. Energy thresholds for p, d, t 
in the CsI(Tl) detectors were 8, 11, 14 MeV respectively, 
and for neutrons in the liquid scintillators 2.0 MeV. De- 
tails about the experimental setup and the particle en- 
ergy determination are given in Refs. [17,19,20]. 

Fig. 1 shows the ratios of the n, p, d and t kinetic 
energy yields measured in 36 Ar+ 124 Sn and 36 Ar+ 112 Sn 
(note the different scale in the figure for n as compared to 
p, d and t). An equal number of events is sorted for the 
two Sn-targets. The different solid angle coverage of n 
and p detectors is accounted for, and the neutron energy 
is efficiency corrected [17]. One can notice not only a sub- 
stantial enhancement of the n yield for the neutron-rich 
system (as may be expected) , but also that the enhance- 
ment is strongly energy dependent. Furthermore, the p 
yield is reduced at low energies for the neutron-rich sys- 
tem, and the t yield is enhanced over the whole energy 
range. On the other hand, the yields of the deuteron 
spectra are the same for the two systems. 

The correlation function, C(q,P tot ) = kN c (q,P to t) / 
N nc (q, Ptot) , is constructed by dividing the coincidence 
yield (N c ) by the yield of uncorrelated events (N nc ) con- 
structed from the product of the singles distributions [18] . 
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q = — Vijm-i), where u is the reduced mass, is 

the relative momentum, and P to t = Pi + Pi is the total 
momentum of the particle pair. The correlation function 
is normalized to unity at large values of q, 80<i7<120 
MeV/c for pp and np, and 160<(7<200 MeV/c for all 
other particle pairs. 
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FIG. 1. Ratios of the n, p, d and t kinetic energy yields 
measured in 36 Ar+ 124 Sn and 36 Ar+ 112 Sn. The ratios are ar- 
bitrarily shifted in the y-axis (by 5 units for n, by one unit 
for p, d, t). The dashed lines correspond to unitary ratios (at 
the angles indicated in the figure). 

The 54°<6><114/120° pp/np correlation functions are 
presented in Figs. 2a, b. The neutron energy threshold is 
here 8 MeV, to match the proton threshold. The shape 
of the correlation functions looks as expected from the 
interplay of quantum statistical and final state interac- 
tions. Comparing the two Sn-targets, one observes a 
small but hardly significant enhancement of the corre- 
lation strength for 124 Sn, in both pp and np correlation 
functions. For the interpretation of the correlation data, 
it is important to note that the correlation function de- 
pends on the space-time extent of the emitting source. 
From the size of the source, a stronger correlation is ex- 
pected for the smaller 36 Ar+ 112 Sn system, an effect ex- 
pected also because of the larger excitation energy per 
particle available for this system (yielding a shorter emis- 
sion time). On the other hand, the change in neutron 
number implies a different symmetry energy which also 



affects the n (and p) emission times. Neutrons are ex- 
pected to be emitted faster in the neutron-rich system, 
which would lead to an enhancement of the correlation 
strength for 36 Ar+ 124 Sn. Thus, the net influence on the 
correlation function is not easily predictable, both due 
to the uncertainty in the symmetry energy and to the 
presence of more than one source of emission. 
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FIG. 2. Angle-integrated (54°<6»<114/120°) pp (a,c,e) and 
np (b,d,f) correlation functions from 36 Ar+ 112 Sn (filled cir- 
cles) and 36 Ar+ 124 Sn (open circles). Low-Pt ot (c,d) and 
high-P t0 i gated (e,f) correlation functions are also shown. 

The emission of light particles from 61 A MeV (mid- 
peripheral) heavy ion reactions originates from (at least) 
three sources, a projectile-like (PLS) and a target-like 
(TLS) evaporative sources (statistical evaporation, SE) 
and an intermediate velocity source (IS). The IS source 
represents dynamical emission (DE), which is described 
either by early nuclcon-nuclcon collisions or by other pre- 
cquilibrium processes. Some evidence for the creation 
of a "neck-source" with low density [21-29] has recently 
been presented, and this may well comprise the domi- 
nating part of the DE source. By selecting nucleon pairs 
with high (low) total momentum, in the proper emission 
source frame, the DE (SE) emission should be favored. In 
the following, the P to t distributions for pp and np pairs 
are calculated in the frame of forward moving sources 
of velocity 0.02c (TLS) or 0.18c (IS). These velocities 
are estimated from Maxwell-Boltzmann fits to the sin- 
gles energy distributions. Low- and high-P tot selections 
correspond to 20-35% of the total correlation yields. 

Figs. 2c, d present the pp and np correlation functions 
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gated on low- Ptot- The pp correlation function is sen- 
sitive to this gate, and a suppression in the correlation 
strength is observed, while the np correlation function is 
only slightly suppressed. This is a first hint that the np 
correlation function may be dominated by DE. 

Figs. 2e,f present the pp and np correlation functions 
gated on high-P tot . Both pp and np correlation functions 
are affected and, in accordance with earlier observations 
[18], an enhancement in the correlation strength is ob- 
served relative to the ungated correlation functions. 
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FIG. 3. pp correlation functions (a,c,e) measured in the 
forward detectors 30° <6><42° , and np correlation functions 
(b,d,f) measured in the range 54°<6><120°, with E n >2 
MeV and E p >8 MeV, from 36 Ar+ 112 Sn (filled circles) and 
36 Ar+ 124 Sn (open circles). Low-Pt ot (c,d) and high-P tot gated 
(e,f) correlation functions are also shown. 

All results from these angle integrated data (54°- 
114/120°) show smaller isospin effects than could be ex- 
pected from the N/Z ratio, which is 1.18 for 36 Ar+ 112 Sn 
and 1.35 for 36 Ar+ 124 Sn. In case the neck formation 
dominates the DE process, this is surprising, since the 
neck region should exhibit a strengthened neutron ex- 
cess, quite sensitive to the isospin of the target nucleus 
[4]. In an attempt to select phase-space regions that fa- 
vor DE emission, we now look at: i) pp correlations in 
the forward (30°<6'<42 o ) region and ii) np correlations 
which include very low energy neutrons (E n >2 MeV, in 
the lab system). Neutrons with 2-8 MeV energy, de- 
tected at 54°<6><120°, correspond to 15-45 MeV energy 
in an IS (u=0.18c) source. 

The results for pp and np correlation functions ob- 



tained under conditions (i) and (ii), are shown in Fig. 3. 
A stronger pp correlation function is now observed for the 
more neutron-rich system, 36 Ar+ Sn. This is true for 
integrated (a) as well as low-P tot gated (c) correlations. 
The low-P tot gated correlations show the largest isospin 
signal. The reason for the relatively large value of the Ptot 
cut (450 MeV/c), is the necessity to populate the large q 
(normalizing) region in the forward (30°<6><42°) detec- 
tors. This cut corresponds to TLS protons with P p <30 
MeV, and IS protons with E p <20 MeV. It is these low 
energy protons that show the largest isospin effect. 

The np correlation functions at angles 54°-120° (Figs. 
3b,d,f), show large isospin effects when low energy neu- 
trons, 2<P n <8 MeV, are included (compare with Figs. 
2b,d,f). Again, the correlation is stronger for the more 
neutron-rich system, indicating that the space-time ex- 
tent of the emitting source is smaller in 124 Sn. In partic- 
ular, both low-P tot gated (here P tot < 240 MeV/c) cor- 
relation functions (Fig. 3d) are enhanced as compared 
to Fig. 3b (instead of being suppressed as expected for 
an evaporative source). A plausible explanation for this 
is that the pairs selected by this gate are dominated by 
DE particles. Since this effect did not occur in the data 
where the low energy neutrons were removed (Fig. 2d), 
we attribute the correlation enhancement to DE parti- 
cles in the energy range 15-45 MeV, which are favored 
in the neutron-rich system. If this isospin effect mainly 
comes from the intermediate energy particles, it should 
be connected to the density dependence of the nuclear 
symmetry energy [15,16]. 

Within the multi-source reaction mechanism described 
above, composite particles, like deuterons and tritons, are 
believed to be predominantly emitted from the DE source 
[30,31], where they are formed by a coalescence mecha- 
nism [32] . Let us first remark that neither integrated nor 
Ptot gated dd, tt and dt correlation functions show any 
appreciable difference between the two Sn-targets [33]. 
This is in agreement with the small sensitivity shown 
by the calculations of Ref. [34]. Even so, a variation of 
the correlation functions such as nd, nt, etc. may be ex- 
pected, as a consequence of the isospin effects on neutrons 
and protons. Indeed, this is the case in our experimental 
data, which can be taken as a further evidence for the 
presence of true isospin effects. Figs. 4a, b present the 
angle-integrated (54 o <6*<120°) nd (a) and nt (b) corre- 
lation functions, measured for the two Sn-targets, with a 
neutron energy threshold of 2 MeV. The anticorrelation 
observed for nd pairs (Fig. 4a), has been observed earlier 
for smaller collision systems, both at lower energies [35] 
and for 61 A MeV [19]. It may originate from the de- 
pletion of low relative momentum nd pairs due to triton 
formation [36]. In Fig. 4a the difference in correlation 
strength between the two Sn-targets can be observed. 

The correlation functions of nt pairs (Fig. 4b), exhibit 
a broad peak which contains the contributions from the 
particle-unbound ground state of 4 H, and possibly from 
higher lying excited states [37]. Once again, a small 
enhancement in the strength is observed for the 124 Sn- 
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FIG. 4. From 36 Ar+ 112 Sn (filled symbols) and 36 Ar+ 124 Sn 
(open symbols): (a,b) nd, nt (54°<6K120°); (c,d) pd, pt 
(30°<6><42°). 

The pd and pt correlation functions, measured in the 
forward angular range (30°<6><42°), are shown in Figs. 
4c, d. The pd correlation function (c) is characterized by 
a pronounced anticorrelation at small q, due to final state 
Coulomb repulsion. The isospin effect is negligible. The 
pt correlation function (d), contains resonance contribu- 
tions from several excited states of 4 He [37]. A small 
isospin effect is seen in this correlation function. 

In summary, isospin effects have been investigated in 
the E/A = 61 MeV 36 Ar+ 112 Sn, 124 Sn reactions, and, 
for the first time, correlation functions from systems sim- 
ilar in size but with different isospin have been experi- 
mentally determined. Both angle-integrated and total- 
momentum gated correlation functions for all different 
pairs of particles containing n, p, d and t, have been mea- 
sured. The largest effects from the isospin of the emitting 
system are seen in the np correlation function. In partic- 
ular, gated np correlation functions which should favor 
a dynamical emission source, show a stronger correlation 
for 124 Sn than 112 Sn. This could be explained by different 
time distributions, with a shorter average emission time 



for the neutron-rich system. Smaller isospin effects are 
also seen in pp, pt, nd and nt correlation functions. These 
experimental results demonstrate that two-particle corre- 
lation functions indeed provide an additional observable 
to probe the isospin dependence of the nuclear EOS. 
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